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Abstract. As it is well known, the classical works of the Dark Current Spectroscopy 
method allow - using some not too accurate theoretical relations, but huge numbers of 
dark current values for thousands of pixels - the evaluation of a reduced number of basic 
impurities parameters. Unlike these works, this paper tries to obtain — by means of some 
better approximations of the Shockley-Read-Hall (SRH) model - more information about 
the studied impurities, as well as the study of the compatibility of the used theoretical 
model SRH relative to the experimental data. In this manner, both the compatibility SRH 
model with the studied experimental data was checked up, as well as the values of some 
additional physical parameters of the impurified semiconductor (the logarithms of the 
pre-exponential factors InDiff, InDep, the effective value of the energy gap Eg) and of the 
Separate capture cross-sections Oy, 0), of the free electrons, and holes, respectively, by 


the studied deep-level contaminants, were evaluated. 


Keywords: Dark Current Spectroscopy, Charge Coupled Devices chips, Dark Current, Shockley- 
Read-Hall model, Deep-level impurities, Capture Cross-Sections of Free electrons and holes 


1. Introduction 


As it is known (but not always checked up, by means of some specific 
compatibility criteria), the dark current in CCDs is described by the quantum 
theoretical model of Shockley-Read-Hall (SRH). It results that: 


a) the dark current in CCDs in described by a huge number (unlimited, 
apparently) of (independent) uniqueness parameters, i.e. the studied CCDs are 
COMPLEX SYSTEMS, 


b) this imposes the use of some APPROXIMATE RELATIONS (as the 
Arrhenius' one or the relations used by the classical works of the Dark Current 
Spectroscopy), 


c) the values of the physical parameters obtained by means of these 
approximate relations depend on the specific features of each pixel, hence they are 
EFFECTIVE VALUES corresponding to some EFFECTIVE PARAMETERS, as: 
(i) the activation energy EF, from the Arrhenius' relation, (ii) the effective energy 
gap E, intervening in the approximate relations, etc. 
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d) in such conditions, it is possible to intervene some "hidden" relations, as 
the Meyer-Neldel one, but also (with different values of the correlation 
coefficient) some additional ones relative to certain parameters of the diffusion 
and depletion dark current, respectively, etc., the study of these (hidden) co- 
relations representing one of main tasks of this work. 


In order to study the compatibility of the quantum SRH theoretical model (SRH) 
[1], [2] with the experimental data referring to the temperature dependence of the 
dark current in Charge Coupled Devices (CCD), starting from the experimental 
results reported by us in the frame of the works [3] - [5], we studied the evaluation 
of the corresponding dominant uniqueness parameters [6]. 


We have found that the minimal set of uniqueness parameters which ensure a 
sufficiently accurate description of the temperature dependence of the dark current 
in CCDs corresponds to: a) the logarithms of the pre-exponential factors 


In Deo,dift» 1nDeO,dep Of the diffusion and depletion dark current, respectively, b) 


the energy gap E, of silicon, c) modulus |Et-Ei| of the difference of energies 
corresponding to the capture traps (of free electrons or holes) inside Si, and to the: 
d) so-called “polarization degree” d of the capture cross-sections of free electrons 


On-Op 
o, and holes o, , defined as: d = arg tanh, ———— }. (1) 
On top 


The corresponding HSR expression of the total dark current was written as (see 


[5], [6]): 


E 
De” (T) = Degg (1) + Degep(T)=T* ex Deo diff - * + 


a E E, -E; 
+7 3/2 | ed ae ai A | (2) 


The accomplished study pointed out the compatibility of the HSR model 
with the indicated experimental data and allowed the evaluation of the chosen 
dominant uniqueness parameters, the obtained values being in agreement with the 
existing ones, obtained by means of different experimental methods [7], [8]. 


2. The Charge Coupled Devices (CCDs) as complex systems 


In the frame of the study [5], we have found that even the classical (HSR) description of 
the CCDs semiconductor material requires a huge (unlimited, practically) number of 
uniqueness parameters (of usual — symbols): Dns X¢> Apix.NA,Me,mp, 


Eas Magis Ops OnsVaa Np Pe |E,-E;|, etc, many of them [e.g. n, p, nj, etc, 


being also temperature dependent, hence introducing some additional uniqueness 
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parameters, as y,E.,E,, etc., which are also temperature dependent, implying 
other uniqueness parameters, and so on]. 


Additionally, the electrons transitions from their “condensed” state in the valence 
band towards the free (“gaseous”) state in the conduction band can be seen as a 
phase transition. For this reason, some descriptions of the Arrhenius’ type of the 
temperature dependence of the currents in semiconductors are to be expected. 


Finally, a CCD is composed by a huge number (of the magnitude order of 10°) of 
pixels, with different and randomly distributed physical properties (see e.g. [9] 
and Fig. 1), which imposes a statistical approach of their features, hence all 3 
basic characteristics of complex systems: huge number of uniqueness parameters, 
phase transitions and necessary statistical descriptions are reunited [10]. 
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Fig. 1. Dark current normalized to 100 e7/s [9] 


3. The dominant uniqueness parameters of the temperature dependence of 
electronic currents as effective parameters averaged over temperature 


As it is well-known (see e.g. [11]), besides the directly measurable parameters, 
there is a large category of parameters whose values can be estimated starting 
from certain (assumed as valid) theoretical relations — the so-called effective 
parameters. 


3.1. Thermionic Emission 


Even in the rather simple case of the thermionic emission of metals 
(sometimes covered by a rather thin oxide layer), the description of the 
temperature dependence of the thermionic current requires the use of some 
effective parameters (averaged over temperatures) — the dominant ones being the 
pre-exponential factor (thermionic constant) A and the work function (extraction 
energy) g defined by means of the Richardson-Dushman relation [12]: 
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: 2 Pp 
T)=A:T~ -exp| -— 
js(T) (2). 3) 
where j,(7) is the saturation current density at temperature T and k is the 
Boltzmann’s constant. 


The classical work [13] presents (in the frame of table 4 of Chapter 3) a collection 
(starting from carbon: A =~ 30 A-cm™-K”, g ~ 4.34 eV, up to uranium: 
~ 6 A-cm”’:K”, 3.27 + 0.05 eV) of 23 concomitantly estimated (by means of the 
least-squares fit) pairs of constants A and g, for different metals. The correlation 
coefficient corresponding to these pairs is r ~ 0.2275 (and the square mean 
relative deviation corresponding to this regression line s ~ 529.07%), hence the 
effective uniqueness parameters A and g of the thermionic emission are 
independent. 


3.2. Dark Current in Charge Coupled Devices 


The dominant uniqueness parameters of the temperature dependence of the dark 
current in CCDs are obtained by means of some: a) partial, b) general (total) 
averages over temperatures. Unlike the state parameters InDiff(T), E,(T), etc., 


which do not depend on the impurities features (concentrations, cross-sections, 
etc), the effective parameters In Diff, Eg ef., Eq, etc. depend on these features and 


on the considered pixel, implicitly. 


a) The general (Arrhenius’ type) CCDs dark current parameters averaged 


over temperatures 


Taking into account the above (negative) numerical result concerning the 
correlation of the pre-exponential factor A and of the corresponding work function 
(extraction energy) g of the thermionic emission of metals, it results that almost 
sure the very strong (see e.g. [9] and Table 1 in following) Meyer-Neldel’s type 
correlation [14] between the pre-exponential factor De, and the Arrhenius 
activation energy E, of the CCDs dark current [15]: 


De(T) = De, cx - ie) (4) 


corresponds to a true physical relation between these effective parameters. 
Because this relation is not obvious, we will name in following such co-relations 
as « hidden » ones [16], they corresponding so to complex systems, with a huge 
number of (apparently) uniqueness parameters, sometimes in rather strong 
relations. We will underline here that our study [17] of several types of numerical 
simulations of different physical processes did not point out any type of numerical 
phenomenon (intervening in the least-squares procedures) which could be 
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misleading about the true or apparent (artifact) character of the studied physical 
relations, indicated by the statistical correlation coefficient. 

b) Partial (SRH type) CCDs dark current parameters averaged over 
temperatures 


By means of a similar procedure, we can easily find that all 5 dominant 
uniqueness parameters intervening in the expression (2) of the dark current for the 


HSR model: In Deg gige = In Diff , InDeo dep, = In Dep , Eg eff, =Eg, E,—E; and d are 


effective parameters, whose values depend on the considered pixel, but are 
averaged over temperatures. Particularly, the HSR uniqueness parameter 
Eg eff. =Eg has a net distinct physical meaning in comparison with any of the 


usual energy gap E g (T) parameters. 


4. Study of the (co-)relations between the main uniqueness parameters of 
Dark Current in CCDs 


The main studied uniqueness parameters were: a) the natural logarithms 
corresponding to the pre-exponential factors of the: (1) Arrhenius relation 
describing the temperature dependence of the dark current in CCDs, determined 
by means of the least-squares fit (regression line) method (denoted by InDArrh), 
or as the intercept with the In De axis of the straight-line joining the representative 


points (from the plane {i De, rab corresponding to the extreme (222 and 291 K, 


respectively) temperatures (denoted as InArrh), (11) diffusion (InDiff) and 

depletion (InDep) terms, respectively, of the HSR expression (1) of the dark 

current in CCDs, b) the Arrhenius activation energy determined as the slope of the 
1 


least-squares fit (regression) straight-line: In De = (Z| [EaArrh], or of the 


straight-line joining the representative points at the extreme temperatures [EaLin], 
c) the width of the forbidden band (energy gap) of the studied CCDs pixel (Eg), d) 
the modulus of the difference |Et-Ei| of energies corresponding to the capture 
center (trap) [Ef] and to the intrinsic Fermi level (7), respectively, e) the X- 
coordinate of the studied pixel (Xpixel), f) the Y-coordinate of the studied pixel 
(Ypixel), g) the distance R of the considered pixel from the center of the studied 
CCD chip (Rpixel). 


In order to find also : (i) the independent or co-related character of the considered 
uniqueness parameters, (ii) the influence of the theoretical description choice, we 
have evaluated (see Table 1) the corresponding correlation coefficients for the 
main pairs of the above indicated uniqueness parameters, corresponding to 3 
simplifying hypotheses: a) null values of the modulus |E?t-Ei| and of the capture 
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cross-sections polarization degree d, b) null value of d and considerably larger 


| Et — Ei 


than | value of the ratio | , c) null value of d. 


Table 1. Study of the main co-relations between the dominant uniqueness parameters of the dark 
current in CCDs 


Correlation Simplifying hypotheses 
ereee Et= Ei andd=0 | |Et—Ei|>> kT andd=0 d=0 
In DArrh, EaArrh 0.999918 0.9999656 0.9999022 
In Arrh, EaLin 0.999949 0.9999591 0.9999350 
In DArrh, Eg 0.8106086 0.2403149 0.1353221 
In Arrh, Eg 0.7243172 0.2414372 0.2410846 
In Diff, Eg 0.999945 0.999932 0.9997006 
In Dep, Eg 0.737485 0.8810 0.574449 
In Diff, |Et — Ei - 0.9518593 0.65158 
In Dep, |Et — Ei - 0.9668592 0.747427 
Eg, |Et — Eil - 0.9503131 0.6506456 
Eg, Xpixel - 0.43671 - 0.2813439 0.17089 
Eg, Ypixel - 0.13286 - 0.0976788 0.1269467 
Eg, Rpixel 3.54x1077 1947x1077 5367x107 


The analysis of the numerical results synthesized by Table 1 points out that: a) 
there is not any correlation between the energy width (gap) Eg of the forbidden 
band (as a representative parameter of the physical properties of the 
semiconductor material) and the co-ordinates OX, OY or the distance R to the 
center of the CCD chip, hence the semiconductor properties are randomly 
distributed for the pixels along these axes and around the center O, 


b) the very high values of the correlation coefficient corresponding to the 
pairs (In DArrh, EaArrh) and (1n Arrh, EaLin) indicate that the genuine Meyer-Neldel 


relations [14]: 


(5) 


e ee 
In Deg =\InDeogo+ 
EMN 


associated to the Arrhenius’ relations [15]: 
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= = E 
De =De, exp - e) (6) 


are fulfilled with high accuracy both for the equivalent Arrhenius’ parameters 
(In DArrh, EaArrh) determined as the intercept and the slope of the least-squares fit 
(regression) line: InDe™ = f (=) and as the intercept and slope (In Arrh, EaLin) of 
the straight-line joining the representative points corresponding to the extreme 
temperatures (222 and 291 K), 


c) there is a strong asymmetry between the diffusion and depletion Meyer- 
Neldel’s type relations for the semiconductor energy gap Eg: while the correlation 
(In Diff, Eg ) is very strong, the correlation (In Dep, Eg ) presents a medium or even 
a weak intensity, 


d) the co-relations (In DArrh, Eg ) and (In Arrh, Eg ) keep (at limit) a medium 
intensity only for the obviously inaccurate assumption: Et # Ei, but they disappear 
totally for the calculations corresponding to non-null values of |Et - Eil, 


e) the modulus |Er - Ei] of the energies corresponding to the capture 
centers (traps) Et and to the intrinsic Fermi level Ei, respectively, presents at least 
weak co-relations with all main uniqueness parameters: InDiff, InDep and Eg of 
the dark current in CCDs. 


5. Interpretation possibilities of dark current non-uniformity in CCD chips 
by means of the « hidden » (of the Meyer-Neldel’s type) correlations 


5.1. Meyer-Neldel’s type relations referring to the diffusion dark current 
Taking into account that the correlation (In Diff, Eg ) is much stronger than the 
(in Diff, |Er—Ei]) one [r(In Diff, Eg) ~ 0.999706, while r(in Diff,|Er-zi|) is only 


0.65158], we will neglect the influence of |Et - Ei| values on the diffusion dark 
current non-uniformity. 


The regression line corresponding to the strong (InDiff, Eg) correlation is 
described by the equation: In Diff =i+s-E,, where the intercept (coordinate of the 
crossing point with the InDiff axis) is i = -10.01989 and the slope s ~ 38.3229eV", 
its accuracy being also very high (standard relative deviation of only 0.21048%). 


According to relation (2), the temperature dependence of the diffusion dark 
current can be written as: 


E 
var-afentetronfettiyns 0 
oO 


118 Ralf Widenhorn, Ionel Tunaru, Eric Bodegom, Dan Iordache 


its characteristic temperature being: 7, = - ~ 302.8086 K. 
“S 


Because the characteristic temperature 7, is rather near to the studied temperatures 


(222 ... 291 K) [hence the differences —-= are rather small] and the correlation 


oO 


coefficient r( in Diff, |Et - Ei|) is very high, the non-uniformity of the dark current at high 


temperatures (where the diffusion dark current prevails) is reduced (see Figs. 1). 
5.2. “Hidden” (of Meyer-Neldel’s type) co-relations referring to the 
depletion dark current 


Taking into account the rather near values of the correlation coefficients 
r(In Dep, Eg ) ~ 0.57449 and r(in Dep, |Et—Ei|) ~ 0.747427, we studied the double 


linear regression: In Dep = f(Eg,| Et—Ei|)= i'+s|-Eg+s9-|Et-Eil|, (8) 


determining its basic parameters: 
a) correlation coefficient r( In Dep; Eg,| Et —Ei|) ~ 0.756385, 
b) standard relative deviation ~ 6.01144%, 
c) coordinate of the crossing point with the InDep axis: i’ ~ 0.45434, 
d) slope relative to Ey: s; ~ 13.04339 eV", 
e) slope relative to |E,-Ej|: s7 ~ 49.390135 eV". 


Starting from relation (2), one finds that the expression of the temperature 
dependence of the depletion dark current is: 


I 73/2 0) 
De Jon = exp) i'+E [s ae |He SS 
dep gl] ol UT | t i| 2 IE, -E,| 
cosh} ———— 
kT 

In order to define also a characteristic temperature of the (1nDep,|Er-Eil) 

: : ; : E,-E; 
correlation, we will consider the particular case: oo, when the 


1 


hyperbolic cosine can be approximated as: cosh x= >-e". In this approximation, 


the above relation can be written as: 


E E,-E; 
Deinp #200 F & 1 +. t | 1 +) }e ; (10) 


2k\T T 
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where the characteristic temperatures corresponding to the (In Dep, Eg ) and to the 


(in Dep, |Et- Ei] ) correlation, respectively, are: 
T, =-——- = 444.8426 K, TM = _1-- = 234.956 K. 
2k +s} k-s9 
Given being that: 


a) while the characteristic temperature of the (in Dep, |Et—Ei|) correlation is located 


inside the studied temperature interval (222 ... 291 K), that corresponding to the 
(In Dep, Eg ) correlation is (444.8426 K) rather distant relative to this interval, b) 


the correlation coefficient r(1n Dep, |Et— Ei] ) ~ 0.747427 is considerably larger that of 
the (In Dep, Eg ) correlation: r(In Dep, Eg ) ~ 0.574449, it results that the large 


spreading of the depletion dark current (which prevails at low temperatures, see 
Figs. 1) is due mainly to the weak (In Dep, Eg ) correlation (see also [18]). 


5.3. Implications of the asymmetry of the diffusion and depletion 


Meyer-Neldel’s type relations on the non-uniformity of dark currents 
Assuming that the pairs of individual values (In Dep, Eg ), (1n Dep, |Et—Eil ), 


(In Diff, Eg ), etc. are normally distributed around their average values <InDep>, 


<Eg>, etc. and denoting by X=Eg-<Eg>, Y=InDep—<InDep>, etc. the 
corresponding deviations, the equation of the confidence ellipses is: 


2 2 


x Y x Y 


pea ge, 11 
a(x) or) CX) a)” ~ 


where o(X),o(Y) and r are the corresponding standard deviations and the 
correlation coefficient, while p is a parameter related to the confidence level 
associated to the considered confidence ellipse (11): 


p=-21-r2} maz). (12) 
Scaling the physical units to have equal standard (square mean) deviations: 


o(X), o(Y), it results (see e.g. [11], p. 40) that the ratio of the semi-minor axis b to 


the semi-major axis a of the confidence ellipse relative to its symmetry axes (see 
also Fig. 2) is: 


a (13) 


De. l-r 
a l+r 


One finds so that for the (InDep, Eg ) and (In Diff, Eg ) co-relations, the 
values of this ratio are: 
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(2) = 0.51989 and (2) = 0.012236. 
4 /)\n Dep,Eg 4 /\nDiff,Eg 


These values explain also the considerably higher non-uniformity of the 
dark current at low temperatures (when the depletion dark current prevails) than at 
high temperatures (in conditions of diffusion dark current prevalence). 


InDiff 


Correlation Coefficient 
0.99970 


b 


Spreading|o 
b ~0,012236 
InDiffWalues 
Ey Kg E, 


Fig. 2. The confidence ellipses associated to assumed normally distributed pairs of 
individual values (In Dep, Eg ) and (In Diff, Eg ), for the o(X)=o(Y) scaling. 


6. Conclusions 


The accomplished study of the possible co-relations between the basic uniqueness 
parameters (InDeoqaig¢= InDiff, InDeddep = InDep, E, and |E,-Ej)) of the 
temperature dependence of the dark current in CCDs pointed out the presence of 
some strong or at least weak co-relations relating the (natural logarithms of) pre- 
exponential factors InDiff and InDep both with the energy gap Eg and with the 
modulus of the energy difference corresponding to the capture traps Et and to the 
intrinsic Fermi level Ei. 


This study pointed out also that the experimental finding referring to the 
considerably larger non-uniformity of the dark current at low temperatures than at 
higher ones is due both to: a) the considerably stronger Meyer-Neldel type 
correlation (In Diff, Eg ) than the (In Dep, Eg ) one [see figs. 2], b) the location of 


the characteristic temperatures of the (In Diff, Eg ) and (n Dep, |Et- Ei] ) co-relations 


very near or even inside the studied temperatures interval (which minimizes also 
the corresponding spreading of the dark current values), while the characteristic 
temperature of the (InDep, Eg ) correlation is located considerably outside the 
studied temperatures interval, leading also to a considerably larger dark current 
non-uniformity at low temperatures (when the depletion dark current is prevalent). 
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The obtained results underline so the theoretical and practical importance of the 
study of the “hidden” co-relations between the physical parameters of complex 
materials, being of particular interest for the technical applications of CCDs (see 
[19], [20]). 
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